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ABSTRACT 

The properties of an extension of the New Massive 3D Gravity by scalar matter with Higgs-like 
self-interaction are investigated. Its peturbative unitarity consistency is verified for a family of 
cosmological Bounce solutions found by the superpotential method. They correspond to the lower 
bound A = — 1 of the BHT unitarity window and describe eternally accelerated 3D Universe between 
two initial/final stable dS^ vacua states. 
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The perturbative quantization of the Einstein gravity is known to introduce at one-loop order in 
the gravitational coupling k 2 = \6itG "quadratic" counterterms C c t = 

which cure its UV divergencies but violate unitarity (for d > 4) [1], [2]. In three dimensions the 
problem of perturbative consistency of quantum gravity model including specific "high derivatives" 
terms: 

Snmg(9^ Ao, A) = ^ J dx 3 ^[eR + 64X K - 2A] , K = R^R^ - ^R 2 (1) 

called New Massive Gravity (NMG) was recently reexamined by Bergshoeff, Hohm and Townsend 
(BHT) [3], [4]. It turns out that it is super-renormalizable and unitary (ghost free) under certain 
restrictions on the values of the cosmological constant A and of the new scale (mass) parameter Ao 
for the both choices e = ±1 of the "right" or "wrong" sign of the R-term [4], [5], [6]. The BHT 
- model (1) has two propagating degrees of freedom and unlike 3D Einstein gravity (i.e. Ao = 
case) it admits interesting (vacuum) solutions - black holes, gravitational waves, etc. [7], [8]. Its 
Newtonian limit for e = — 1 reproduces an attractive gravitational force [9]. 

The New Massive Gravity (1) can be realized perturbatively as a massive spin two quantum field 
hav (dfiu = 9fj,u + K>huv) - theory with unusual kinetic term and specific self-interactions in external 
constant curvature background g^ u . Such a procedure provides well defined rules for studying the 
quantum graviton (h^) properties and eventually it might shed some light on the nature of the 
"quantum 3D geometry" behind it. Thus we have an example of consistent (3D) quantum gravity 
model that shares many of the (desired) properties of AD gravity. It is therefore natural to address 
the question about the perturbative quantum consistency of the New Massive Gravity coupled to 
matter fields. Since we are interested to study the Cosmological consequences of such BHT-matter 
model, the simplest and most representative candidate of matter QFT is the one of scalar field 
a(x^) with Higgs - like self-interaction: 

C m att(<T,g^h,ma) = V=g(-\d^a-V(a)) , V(a) = J(a 2 - - 2A (2) 

The parameters 7, m a and A above denote the renormalized coupling constant 7, "Higgs mass" ra a 
and the vacuum energy A and they are assumed to depend on the "graviton mass" scale Ao = ■ 
Since in four dimensions (in curved space) this model is known to be renormalizable [10], [11] one 
expects that it is power - counting super-renormalizable in three dimensions. Note that the relevant 
(lowest order in 7) counter-terms to £ m att are again in the form aR 2 + fiR^R^ [10], i.e. they do 
coincide with the gravitational ones (in the case of minimal coupling £ = 0). We further assume 
that the corresponding (vacuum) "semiclassical" equations [10], [12]: 

Rnv - ^9^R + A<? M * = ^k 2 < vac\f^(h^,a)\vac > (3) 

for this BHT - matter model (to lowest order in the both couplings k 2 and 7) with the UV counter- 
terms added and when adiabatic vacuum states \vac >= \g^ c ,cr± >, (V'(a±) = 0) exist, are 
represented by the classical equations obtained from the following "effective" Matter-Gravity ac- 
tion: 



qeff 
NMGM 



I dx 3 ^g~\eR + 64A A" - ^| Vcr| 2 - V{a)\ (4) 

K J L Z 



,2 



The BHT - cosmological constant A is now included in the potential (2), i.e. V(<j± = = — 2A 
for the case eAo < 0. In order to describe in a compact form the corrections (i.e. backreaction) 
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to the classical background g^v {R^u = Rfiuid^v) etc.) due to the quantum fluctuations of both 
h^ u and a, as usually (see ref. [12]) all the higher derivatives contributions are included in the 
semiclassical "source" in the r.h.s. of eqs.(3). 

The present paper is devoted to the investigation of the properties of homogeneous and isotropic 
spatially flat k = (—,+,+) 3D Universe : 

dsl = -dt 2 + a 2 eff {t)(dx 2 + dy 2 ) (5) 

described by the matter extension (4) of BHT-model (1). Our problem is therefore to find analytic 
solutions of the corresponding 3D "Massive Cosmological" equations 3 : 

a + &ip + V'(a) = 
(p(l + 8eA o 2 ) + ^ 2 (l + 4eA o 2 ) + e(^ 2 -V>)) =0 
e(& 2 + 2V(a))-p> 2 (l + 4e\ <p 2 ) = (6) 

derived from the action (4) above. We are further interested in studying specific features of these 
solutions as for example: the effects of the matter on the BHT - unitarity conditions (u.c.'s), the 
conditions of absence of initial/final singularities, how many different periods of acceleration and/or 
deceleration they describe etc. 

Observe that eqs. (6) for the scale factor a(t) = e 2 and for the scalar field a(t) are of second 
order although the action (4) and the eqs. (3) include (up to) fourth order time derivitives of g^ u . 
The special BHT unitarity motivated choice of the relative coefficients /3 = — |a = 64Ao (that 
fix the form of the K-term for D = 3) together with the FRW anzatz (5) are responsable for the 
cancellation of all the terms involving ip^ 111 ^ and ip( IV ") derivatives and they cure another disease of 
all "higher derivatives" gravities. The family of solutions we are considering does have well defined 
initial values Cauchy problem (i.e. "causality") for all the values of the parameters e, Ao, 7, m 2 , i.e. 
independently of whether or not the BHT u.c.'s [4] are satisfied. The construction of all the solutions 
of this second order system of nonlinear eqs. (6) is however rather difficult problem, although 
numerical methods are indeed available and largely used for this type of dynamical systems. One 
has to be able to further analyze how their properties depend on the initial conditions (do, &o, <p>o, 
(po) and on the values of the parameters e, Ao, 7, vn a as well. 

At this place we take a different route known as "superpotential" method, which provides 
very special analytic (non-perturbative in 7 and Ao) solutions of eqs. (6). It is widely used in 
the construction of domain walls (DW) solutions of different (super) gravity models in arbitrary 
dimensions (see for example refs. [13], [14] for Ao = case). It consists in the introduction of an 
auxiliary function W(a) (called superpotential) such that: 

V(a) = 2eW 2 (l + We\ W 2 ) -2(W') 2 {l + 32e\ W 2 ) 2 
ip = -2eW(a) 

a = 2W'(a)(l + 32e\ W 2 ) (7) 

For Ao 7^ it represents an adapted version of the Low-Zee superpotential [15] introduced in the 
context of DW's solutions of d = 5 Gauss-Bonnet improved gravity. The statement is that for each 
function W(a) that reproduces our matter potential (2) the solutions of the first order system (7) 
are solutions of the eqs. (6) as well. 

3 in what follows we use the convention c = k = Ti = 1 
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We next consider a particular linear superpotential W(a) = Ba , which leads to a family of cr 4 
- like {Z<i symmetric a — > —a) potentials (2) with 

7 = 128A -B 4 (1 - 64A -B 2 ) , ml = -8eB 2 (l - 64A £ 2 ) 

°- - ~W A " ek ~ m2 ■ fnl = V " m ~ -\ V " M - f (8) 

Hence depending on the eAo sign (and whether B 2 < m 2 or B 2 > m 2 ) we have normal or inverted 
double-well potentials for eAo < or the standard a 4 - potentials for eAo > 0. The fact that the 
form of the matter potentials (encoded here in 7 = 7(Ao) and m 2 = m 2 (Ao)) does depend on the 
gravitational mass scale m 2 = is an intrinsic property of the superpotential method (see eqs. 
(7)) when applied for models which actions contain K - type terms (1). 

Note that for the particular choice of W(a) we consider, the parameters {7, m 2 ,m 2 } are not 
independent. As one can see from eqs. (8), they describe a special surface 

2 

ml = 32 7 TO 2 (m 2 + 4e 7 ) , 1 + e|^ > (9) 

in the parameter space of the BHT - matter model (4) . Taking the effective masses m 2 and m 2 as 
independent variables we realize that the coupling constant 7 (or equivalently the "superpotential" 
parameter B 2 ) are given by: 



m 2 /. ml \ m 2 B 2 



7± 




i.e. they take two different values j± (and B±) when < and only one 7_ (or B 2 ^) in the case 

^\ > 0. Together with the signs of m 2 and m 2 they characterize the particular shapes of the 
potential V(a\ r y±,m 2 r ) for which we are able to solve analytically eqs. (6) by the superpotential 
method. We will consider mainly the case fh 2 ^/ < (i.e. eAo < 0) which describes: 

(a) e = — 1, Ao > BHT - model with double- well (Higgs) potential V-(o). 

(b) e = 1, Ao < BHT - model with inverted double well potential V+(a). 

In both cases we have a\ > and therefore o± = ±^/^- are the minima (e = —1) or the 
maxima (e = 1) of the potentials V±. For these "vacuum" values a± of the field a(t), the BHT - 
matter model (4) reduces to the pure BHT - gravity (1) for a particular value of the cosmological 
constant 4 A = m 2 = i.e. we have A < for Ao < (e = 1) and A > for Ao > (e = — 1). 
The effective cosmological constant A e jj is defined by the values of 3D curvature on the 
corresponding vacuum solutions (<7^ c ' =l: , <r±): 

R^Kg^v, <r±) = 2y3 + ^ 2 = -8e(W') 2 (l + 32X eW 2 ) + 6W 2 (11) 

which gives R^ 3 ' \g v ^f ,cr±) = — 16 3 Xq = 6A^. Therefore in both cases when eAo < 0, the A^ = 

2\m 2 \ > is positive, and the scale factor a 2 ac (t) = e 2 V Ae ff l = a 2 dSi represents dSs geometry in 
agreement with the generic relation between A e jj and A [4]: 

Af ff = -2m 2 (e ± y/l + X BHT ) (12) 

We have two options for the initial conditions of a(t) and a(t): 

4 In the BHT notations [4] we have Xbht = — ^ and therefore both eA < cases correspond to the limiting 
value Xbht = — 1. 
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type (1): a 2 (t) < a'±, i.e. 



a(t — > ±00) = cr±, cr_ < a(t) < a + 
at ac = a dS3 (t) = eV 7 ^ t 



type (2) : a 2 (t) > a 2 ± 
- (2a) 



" (2b) 



er_ < a(t) < a(t = to) = 00 
a vac = a dS 3 but a(to) = and R^ 3 \t = to) = 00 

er + < <r(i) < <r(i = to) = 00 



Let us consider first the non — singular solutions corresponding to the initial conditions of the type 
(1) above, i.e. the ones interpolating between the two vacuum states: \vac >= \a±,dSs(A e ff) >. 
We find by direct integration of eqs. (7) two exact solutions: 



<7±(t|e) = 4/^tanh(6 ± (*-to)) , 4(t|e) = e v±(t|e) = ( cosh6 ± (t - t )) S± (13) 



2 7± 

where 5± and 6± are given by: 



b ± = ^=(l±Jl + e^) , 8 ± = -e^ = -2e^ (14) 

Observe that a±(t) formally coincide with the flat space QM "instanton" solution 5 for the double 
well potential V±(a), but now 7 and m 2 are related by eq. (9). The scale factor a±(t) resembles 
the well known bounce solution for closed k = 1 dS% Universe (Ao = 0, a = const) as one can see 
by comparing the behavior of the corresponding Hubble parameters: 

Hi =0 (t,e,X ) = — = -eB±a±(t) , H k= \t) ~ tanh2v / A(t - to) 
a± 

The explicit form of the curvature scalar (11), i.e. 

fl±(t) = cosh^S - £ + k S '" h2 M < - *">) <15> 

confirms that the bounce solution (13) represents non-singular asymptotically dS^ ((a)dSs) geom- 
etry for all t £ (—00, 00), in particular 

R^(t -> ±00) = 12|m 2 | = 6A± / = L 9r » l pl (16) 

Another important feature of this (a)dS^ Bounce solution (13) is encoded in the behavior of decel- 
eration parameter: 



jP_ = -\ s ±( £ ) 

V 2 " smh 2 (b±(t-t )) 



?±(*) --i-2-^ = -i- . — (17) 



3 note that in our curved space example the time is not euclidean (see eq. (5)) so that it just looks like an instanton 
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where s± are real numbers 



Bl \;_ .G J 



2 m 2/ m 2 7 2 V y 2m 2 ^ |m 2 | 

They are determined by the ratio of the two scales 



9r |A e// | 2|m 2 | ' u |m 2 | 

Ipl ^ L gr < L a (18) 

introduced by the graviton mass M 2 r = 2|m 2 | and by the scalar field mass 6 |m 2 |. 

Note that our Bouncing 3D Universe is almost invariant under time reflections T : t — > 2to — t 
(t — > —t for to = 0). The exact invariance is achieved by combining T with the Z2 - charge 
conjugation of a(t): 

C : C<j( fe )(t) = -<T(fc)(t) = (7( afe )(i) 

which interchanges "kink" <J^(t) with the "anti-kink" o~( a fc)(t) = —a^(t) = a^(—t) (i.e. — > cr_ 
in the vacuum sector) . Hence the past of the "kink" coincides with the future of the " anti-kink" . 
This fact indeed reflects the CPT - invariance (Px l = —x l , i = 1,2) of the BHT - matter model 
(4). 

Although the forms of the Bounce solutions (13) in the two cases {e = — l,Ao > 0} and 
{e = 1, Ao < 0} considered above are almost identical 7 , the features of the Universe evolutions they 
describe are completely different. The origin of the qualitatively different physics behind is in the 
form of the matter interaction V±: stable vacuum double well potential for e = — 1 versus unstable 
"maxima" of the inverse double well potential for e = 1 case 8 . Consider first the {e = —1, Ao > 0} 
case. Since s±(e = — 1) > is always positive, the deceleration (17) q±(t\e = — 1) < is negative 
for all t £ (—00,00) and therefore we have an eternally accelerated Universe. It starts its evolution 
from well defined (stable) 0IS3 vacuum state (<r_, (fS^Agjj)) at the left minima of the matter 
potential living a period of (accelerated) expansion (ip > for B < 0) for all t < to. At the moment 
t = to (say to = 0) the Universe enters in the contraction (tp < 0) period towards the new ("final") 
vacuum state (er+, dSs(A^^)) with Act = a + — cr_ = rn a ^J^ and A^ = = 2m 2 . 

In the second case {e = 1, Ao < 0}, that corresponds to the inverted double well potential V-(a), 

B 2 

due to the negative values of s_(e = 1) = — t^t , the deceleration q~(t\e = 1, Ao < 0) changes twice 
its sign. Therefore we have in this case two acceleration epochs (initial and final ones) and one 
period of deceleration. By studying the properties of scale factor a 2 (t) and its derivatives a and a 
we can reconstruct the history of this non-singular Bounce - like 3D Universe. At the infinite past 
t — > —00 when: a(t) — > cr_, a?_(t) — > e 2 !" 1 ^* it starts from the unstable maximum (<j_, ciS^A^)) 
of U_(ct) its first (semi-infinite) period of acceleration and contraction with R^(t\e = l,Ao < 0) 
decreasing from 6A e jj = 2\m I to i? (3) (ti) = 0. At the moment h < 

«i = 




6 or by fh 2 = —^f for the case e — 1, since in this case m 2 < 0, i.e. it describes tachyons. 
7 the only difference is in the scale factors a 2 (e = —1) = a~ 2 (e = 1). 

8 intrinsically related to the wrong/right signs e = ±1 of the Einstein-Hilbert term eR in the action (4) 
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the curvature R^ 3 \ti < t < —t±) < changes its sign entering in the epoch of negative curvature 
of finite duration T_ = 2\t\\ up to the moment i\ = —t\ > when it changes its sign once more 
towards the finite unstable {a + ,dSs} - state. Concerning the accelaration we observe that at the 
moment \t a \ > \t\\: 



t a \ = j- hij^^+yVl + l 



6_ 

the Universe enters in deceleration period which takes place over a finite period of time T^ ec = 2\t a \. 
The second (final) epoch of acceleration and expansion starts at t = \t a \ towards its "future eternity" 
{er + , dSs} at t — > oo. 

An equivalent (but more compact) description of the Universe evolution is based on the so 
called equations of state: p e ff = Peff(peff), where the effective pressure p e ff and energy density 
p e ff are defined by the diagonal r.h.s. of eqs. (3): 

< vac\Tjf(a,h a/3 )\vac >= {~Pef f , Peff , Peff) 

and it appears to be a simple consequence of the eqs. (6), (7): 

Peff(t) = 8em 2 s± + (s± - l)p e ff(t), Pe ff > (19) 

In the {e = — 1, Ao > 0} case we have 

Peff(t) = -8B 2 ± -p eff (t){l-^) , B 2 ± <m 2 (20) 

i.e. p e ff < is always negative and that is why the 3D Universe is eternally accelerated. In the 
second case {e = l,Ao < 0} the pressure changes twice its sign (at the moments ±£ a ), due to the 
existence of critical value of the energy density 

or , 1 x c\\ 8B 2 \m 2 \ 
pf // (6 = l,A <0)= |m2| + jg 2 

where p e ff vanishes. This fact leads to the observation that there exist two periods of negative 
pressure and another one of positive pressure between them. 

An important problem in the classical (and quantum) Cosmology concerns the initial and/or 
final conditions (i/f. c.'s). Namely, whether and how the particular features of the Universe 
evolution depend on the corresponding i/f. c.'s. The BHT - matter model (4) provides an example 
of two distinct families of such conditions (type (1) and (2) above) that lead to qualitatively different 
Universe histories: non-singular Big-Bounce for the type (1) (a 2 < a±) and the singular Big - Bang 
for the type (2) (a 2 > a±) for the same values of all the parameters {k 2 , Ao, m 2 } of the model. The 
corresponding type (2) solutions: 



a(t|e) = y^coth(& ± (t-t )) , a 2 ± (t\e) = \ smh b±(t - t )\ 5± (21) 

are "singular" at t = to, i.e. 

a(t ->• t ) ->• oo, <p(t ->• t ) -> oo, R( 3 \t ->• t ) ->• oo (22) 

and they exist in two disconnected periods of time. For t 6 (— oo,to) we have the same vacuum 
dSs(A e f f)} (stable for e = —1; unstable for e = 1) as starting point. In the case e = —1, Aq > 
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these solutions describe accelerating and contracting to Big-Crunch 3D Universe, while for the 
unstable inverse double-well potential V-, i.e. in {e = l,Ao < 0} case, the initial acceleration 
period ends at the moment t a (t a < to), when q(t) changes sign and the new deceleration epoch 
begins. Similarly the eqs. (21) when considered in the interval to < t < oo represent accelerating 
and expanding Universe created in the moment to from the Big-Bang singularity (R^(to) — > oo) 
towards the infinite future (t — > oo) dSs - vacuum {cr+, dSs(A e f f)} in the "stable" {e = —1, Ao < 0} 
case. Again the {e = 1, Ao < 0} - Universe exhibits periods of acceleration and deceleration. 

The purpose of this short discussion about the effects of the choice of the i/f c.'s {(c?,/, &i,f)} 
and {(<Pi,f, 'Pij)} on the properties of the solutions of eqs. (3) for the BHT - Matter model (4) 
is to demonstrate the incompatibility of the type (2) i/f c.'s, involving scales L2 < l p i (due to 
< >— > 00, < <p >— > 00), with the range of validity L2 ~> l p i of eqs. (3). On the other hand, 
the type (1) i/f. c.'s (i.e. a 2 < a±, (p 2 << 4B 2 a±, etc.) provides an example of "semiclassical 
consistency", when the parameters of the model {e = ±1, A e ff = M 2 r = = j^} are chosen 

such that 9 

(remember L gr = Lry\ is the smallest scale in this approximation). 

The non-perturbative cosmological solutions of BHT - massive gravity coupled to scalar QFT's 
with a 4 self-interaction (4) we have constructed within the restrictions (9), (10) imposed by the 
super-potential method, do not allow a complete discussion of the "unitarity consistency" of these 
models. However one can easily verify whether (and for which values of the parameters e = ±1, Ao 
and m 2 ) the above mentioned restrictions are compatible with the pure BHT - gravity unitarity 
conditions [4]: 

m 2 (A e// - 2em 2 ) > 
Mgr = -em 2 + l -A ef} > A eff (23) 

together with the well known conditions for unitarity and tachyon-free "light" massive scalar cr 4 - 
QFT (2) in dSs classical background : 

< ml < A eff (24) 

We have already realized that the Bounce solutions (13) in both cases eAo < 0, i.e. of normal and 
inverted double well potentials V±(a), correspond to the limiting value Xbht = — ^ = — 1 of the 
BHT "cosmological" parameter Xbht- As a consequence the effective cosmological constant 



A tff = A eff = 21"' 
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is equal to the "graviton (mass) 2 " , i.e. M 2 r = A e ff = 2\m 2 \. In the case: e = —1, Ao > the BHT 



- u.c.'s (23) are automatically satisfied and the restriction (9) now reads 

47 < m 2 a < 2m 2 (25) 

Therefore the BHT - matter model (4) with "wrong" sign e = — 1 of the curvature term, Ao > and 
V + (a) representing double well potential (2) with rhlj < and 7 > is unitary and tachyons-free. 
It corresponds to the lower bound Xbht = — 1 of the dS% unitary window [4]: 

-1 < Xbht < (26) 



9 in both cases the treatement of (sub) Planckian scales L < l p i requires an introduction of other degrees of freedom 
(than h^v) and other methods of quantization 
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This special point of the pure BHT - gravity model (1) is known to describe "partially massless" 
graviton of one propagating mode, instead of the two helicities ± 2 graviton massive modes for the 
generic M 2 r > A e ff case. The observed reduction of the physical degrees of freedom of h^ v (for 
^bht = —1) turns out to be a result of the enhancement of the gauge symmetries of the linearized 
(effectively Pauli-Fierz) form of the corresponding equations of motion [16], [4]. 

We have indeed similar non-singular solutions (13) for the case e = 1, Ao < (i.e. m? < 0) 
representing the inverted double well potential V_(a) (7 < 0, fh 2 > 0). However the BHT - matter 
model for this range of the parameters is neither unitary nor tachyons-free. For e = 1, m 2 < and 
A e jj = 2\m 2 \ the first of the BHT - u.c.'s (23) does not hold and now a± - dS% "vacuum" represent 
the maxima of V-(a). Hence they are unstable and the corresponding asymptotic (t — > ±00) a - 
particle states are tachyonic (m 2 ± < 0). Due to the "right" sign e = 1 of the R - term one can 
try to interprete this model as an effective action for the scalar QFT with V- (a) self-interaction in 
classical gravitational background, thus relaxing the BHT - quantum gravity unitarity conditions 
(23). In order to avoid tachyonic instabilities one should consider "semiclassical" states build on 
the metastable minima of V-(a) (i.e. ao = 0) and to further study their decay. Note that such 
state is qualitatively different from the "constant curvature" \a±, dSs(A e ff) > - states. As one can 
easily see from eqs. (6), now we have 

<r(0) = 2B, 0(0) = -4B 2 , ip(0) = 

As a result this relative minima state describes linearly growing matter and (non-constant curva- 
ture) Gaussian scalar factor 

og(i) = e- 2eBH \ a (t) = 2Bt (27) 

which makes the analytic study of the spectrum of the fluctuations around it rather complicated. 
Let us briefly discuss the remaining two models when eAo > (or equivalently rh 2, y > 0), i.e. 

V, rM = ± Mp + «) 2 -2A (28) 

n\m\ 



The solutions of eqs. (7) are now periodic of period T = B IJ^ 



a cyc (t)= l -^-t & n(-^-t) , a 2 cyc (t)= cos (-j^pt) * (29) 

and are defined on the finite interval t € ( — -j, -j) only, where they represent initial (Big-Bang) and 
final (Big-Crunch) curvature singularities with a cyc (±^) — > ±00 at both "ends". The main problem 
with these models (and of their solutions (29)) is that they do not admit "vacuum" solutions of 
constant curvature. The absolute minima/maxima (for e = ±1) of V± c (o~), do = corresponds to 
singular non-constant curvature solution (27). Note that the BHT - analysis of the perturbative 
consistency of 3D gravity is based on the linearization g^ v = g™ c + nh^ around constant curvature 
stable vacuum states \cr±,gj^ c > (i.e. the ones corresponding to M3, AdSs, or dSs). They are 
indeed present in the pure BHT - gravity (1) for all values of the parameters e, Ao and A. When 
matter is added, the presence or absence of such type of vacuum solutions depend on the form of the 
potential V(a) and on the values of the parameters e, Ao, m 2 . As we have seen on the examples of 
the double-well potentials (eAo < 0) the analysis of perturbative "semiclassical consistency" of the 
BHT - matter models follows the lines of the BHT - arguments [4], [5] with special attention to the 
unitarity and renormalizability of the matter sector. In the eAo > type of models characterized 
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by the potential (28) such an analysis requires new methods and their consistency remains an open 
problem. 

We have avoided up to now the discussion of the important question about the possible manners 
to couple Matter QFT's to the New Massive Gravity To begin with the natural questions to ask 
are: Should one also include higher derivatives terms like : i?|V<r| 2 , |Vo"| 4 , R^ u \7 ^aX7 u a , RV(a) 
etc.? How such terms might influence the unitarity properties of the corresponding Matter-Gravity 
model?, etc. Our particular choice (4) was motivated by the renormalizability of the matter QFT 
and we have included in the effective action the UV counter-terms (all of type K for the minimal 
coupling £ = 0) that appear in the lowest order in the couplings k 2 and 7. The analysis of the 
structure of divergent diagrams of this graviton -scalar QFT (in constant curvature background) at 
orders k 3 and k 4 suggests that new counter-terms like those mentioned above could be present. This 
line of arguments however unavoidably leads to the introduction of large number of new parameters 
in the corresponding effective action. On the other hand, the original unitarity proof of Bergshoeff, 
Hohm and Townsend [3] that is based on the equivalence of the NMG action (1) to the specific 
"second order" derivatives (including one more auxiliary spin two massive field f^ v ) action 

1 r ™2 

eff 



S PF 



= JLJ dx 3 V=g{ -R + rc, v - ^(ru„ - f 2 )} (30) 



(with / = g^ u f^ u ) appoints an economic alternative without any new parameters involved. Namely 
to consider the following rather obvious generalization of the Pauli - Fierz inspired action (30) 
including matter fields: 

Spfm = ^ / dx 3 V=9{eR - ^|Va| 2 + V{a)) + f^E, u - ^V<%, - f 2 )} (31) 
where we denote 

K 2 1 

E ilv = G iXV - —T^ , T IIU {a)=V tl aV u a--g flu V p aV f) a-g llu V(a) (32) 

The matter potential V(a) has to be chosen by renormalizability arguments. Similarly to the pure 
gravity case [3], [17] one can easily eliminate f^ v in order to get the following suggestive form of 
the New Massive Gravity - Matter action: 



S$Lqm = ^J d^V=g{eR + ±K - k 2 (\\V<j\ 2 + V(a) 



k R ^-\ RT ") + ^ 

Observe that the last two (new) terms: 



U = R^T^ -hiT^R^V^VvG-hlV -^|Va| 2 

Y = T^T, V -\{T£) 2 = h -\Vo\* - \V 2 - \V\V0\ 2 (34) 

are of order k 3 and k 4 correspondingly since the k -expansion of both R^ v and R starts by order n 
terms. Hence the above action (33) when terms up to order k 2 only are taken into account reduces 
to the action (4) we have studied. Although the symmetry principles that single out such higher 
derivatives Gravity-Matter actions among the many others (with much more free parameters) are 
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not completely clear their unitarity and renormalizability indeed represent an interesting research 
problem, which however is out of the scopes of the present paper. Let us just mention that the 
unitarity of the matter sector in the complete action (33) is rather problematic due to the presence 
of the both |Vcr| 2 and |Vcr| 4 terms. The propagator of such massive scalar field contains ghost 
contributions unless certain fine tunning of the coefficients takes place. 

The main objectives of the investigations presented in this paper are the properties and the 
perturbative consistency of the BHT - matter models (4) with specific renormalizable scalar 3D 
QFT added. We have considered as an example the a 4 interaction potential that admits degenerate 
adiabatic vacuum states. Our goal was the construction of the explicit non-vacuum solutions of 
the corresponding "effective" action's equations that possess "maximal" 2D Poincare symmetry 
0(2) <8>T2 and that interpolate between two such vacua. The particular Bounce solutions of the 3D 
Massive BHT - Cosmology with scalar matter (2) we have found by superpotential method turns 
out to saturate the lower bound Xbht = — 1 (i-e- h-eff = 2|m 2 |) of the BHT unitarity window. 
There exist strong indications that these solutions are stable as well, due to the fact that they 
carry Z2 topological charges ±1 . The power of the superpotential method we have used in this 
paper allows the construction of analytic solutions for a vast variety of potentials. The simplest 
case of linear superpotential we have studied in detail represents just one example selected by the 
renormalizability of a 4 interactions in three dimensions. The generalization of our results to the 
case of simplest periodic superpotential W(a) = A + Bcosa that corresponds to the modified 
Double Sine-Gordon potential V(o~) is straightforward. Let us also mention the special case of 
BHT- gravity coupled to free massless scalar field , i.e. V(a) = const = — 2A. In this simple model 
a family of unitarity consistent Bounce solutions can be easily found for arbitrary values of the 
cosmological constant < |A| < m 2 without the use of any superpotentials. 

Another important line of investigations concerns the construction of domain wall solutions of 
the Massive 3D Gravity -Matter models (2) and their further applications for the reconstruction 
of the renormalization group flows in the perturbed 2D CFT dual to these BHT-matter models. 
The well known relation between the FRW cosmological solutions and the static DW's permits 
an extention of our results to the case of Janus kink-like DW's for the model (4) that interpolate 
between two stable AdSs vacua. It turns out that they correspond again to the special value 
Xbht = — 1, but now with A. e fj = — 2\m 2 \ [18]. 

In conclusion: the unitarity consistency of the simplest Massive 3D Gravity - Matter model 
(4) we have demonstrated on the particular class of cosmological Bounce solutions confirms the 
expectations that the BHT - model indeed admits physically interesting renormalizable matter 
extensions. 
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